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We study the phase diagram of the continuum traffic flow model of a highway with an on-ramp. 
Using an open boundary condition, traffic states and metastabilities are investigated numerically 
for several representative values of the upstream boundary flux / up and for the whole range of 
the on-ramp flux f mp . An inhomogeneous but time-independent traffic state (standing localized 
cluster state) is found and related to a recently measured traffic state. Due to the density gradient 
near the on-ramp, a novel traffic jam can occur even when the downstream density is below the 
critical density of the usual traffic jam formation in homogeneous highways, and its structure varies 
qualitatively with / rmp . The free flow, the recurring hump (RH) state, and the traffic jam can all 
coexist in a certain metastable region where the free flow can undergo phase transitions either to the 
RH state or to the traffic jam state. We also find two nontrivial analytic solutions. These solutions 
correspond to the standing localized cluster state and the homogeneous congested traffic state (one 
form of the novel traffic jam), which are observed in numerical simulations. 
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I. INTRODUCTION 

Traffic flow, a many body system of strongly interact- 
ing vehicles, shows various complex behaviors. Numer- 
ous empirical data of the highway traffic have been ob- 
tained , which demonstrate the existence of distinct 
dynamic states and dynamic phase transitions between 
them. Recent studies reveal physical phenomena such 
as hysteresis, self-organized criticality, and phase transi- 
tions in the traffic flow |6|^]. 

The transition from the homogeneous free flow to 
the jammed state has been studied by microscopic and 
macroscopic models without any inhomogeneity in the 
system Ml. The traffic jam, one of the dynamic phases 
of the traffic flow, appears spontaneously when the ve- 
hicle density is between the two critical values p c \ and 
Pc2 (> Pel)- The traffic jam, however, can appear even 
below p c i. The traffic jam can be triggered by localized 
perturbations provided that the density is larger than a 
different critical value ph (< Pel)- As a result, in the 
density range between p^ and p c \, both the free flow and 
the traffic jam can exist, resulting in metastability and 
hysteresis [p^|-p7[. It is observed that some features of 
the traffic jam are uniquely determined by underlying 
dynamics, and independent of initial conditions of the 
traffic flow that lead to the jam Q. The presence of such 
characteristic features is also reproduced by analytic and 
numerical studies of traffic flow models [l5|JT^ |. 

The synchronized traffic flow, another dynamic phase 
of the traffic flow, is identified in recent measurements on 
highways The synchronized traffic flow resembles 

the traffic jam in the sense that both states produce in- 
homogeneous density and flow profiles. The dynamics of 
the synchronized flow is however much more complicated 
than that of the traffic jam. One notable property of the 
synchronized traffic flow is the high level of its average 
flow, which almost matches the flow of the free flow state. 



The synchronized traffic flow is observed, in nearly all oc- 
casions, localized near ramps and it is thus believed that 
ramps are important for the stability of the synchronized 
traffic flow. The discontinuous transition from the free 
flow to the synchronized flow can be induced by local- 
ized perturbations of finite amplitudes. Measurements 
show a hysteresis effect in the phase transitions between 
the free flow and the synchronized flow: the transition 
from the synchronized flow to the free flow occurs at a 
lower on-ramp flux, or lower upstream flux, than that for 
the reverse transition. In Ref. ]l9| , the recurring hump 
(RH) state is proposed as an origin of the (nonstation- 
ary type) synchronized traffic flow H, and the dynamic 
phase transitions between the RH state and the free flow 
are investigated using continuum traffic equations that 
take into account the effect of ramps. In the RH state, 
the vehicle density and the velocity show temporal os- 
cillations which are localized near on-ramps. That the 
synchronized flow is maintained for several hours can be 
explained from one important property of the RH state, 
its being a limit cycle of the traffic equations. The RH 
state can be characterized as a self-excited oscillator, 
where constant vehicle flux from an on-ramp serves as 
a source of the repeated excitation and each excitation 
is subsequently relaxed within a localized region. The 
traffic equations also describe the hysteresis phenomena 
between the RH state and the free flow. 

The traffic jam and the synchronized traffic flow are 
distinct phases of traffic flow. However, the distinction 
between the conditions for the appearance of the jam 
state and the synchronized flow is not clearly identified 
yet, both in measurements and in model studies. High- 
way measurements analysis reports that almost identical 
initial states of the traffic flow can evolve to both the 
traffic jam and the synchronized flow ||. 

To describe the hysteretic phase transitions between 
the free flow and the synchronized flow, a different macro- 
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scopic model based on a gas-kinetic approach is also pro- 
posed J2(|. In this model, a peak of the inflow from an 
on-ramp produces a congested but homogeneous region 
near the on-ramp, which spreads in the upstream direc- 
tion. This homogeneous congested traffic (HCT) state 
is proposed as an explanation for the (stationary type) 
synchronized traffic flow ||. The subsequent study jgjj ] 
investigated the phase diagram of the model and identi- 
fied additional dynamic phases such as the standing lo- 
calized cluster (SLC), the triggered stop and go (TSG), 
and the oscillating congested traffic (OCT) states. An- 
alytical conditions for the existence of these phases are 
provided and it is suggested that the phase diagram is 
universal for a class of traffic models. The study is how- 
ever restricted to the traffic states generated from a par- 
ticular initial condition and thus important issues such 
as multistability and hysteresis are not addressed. 

In this paper, we investigate the phase diagram of the 
traffic flow in the presence of an on-ramp using a different 
continuum model [flif , which tests the idea of the univer- 
sal phase diagram. Various traffic states in Ref. [|T| are 
reproduced. However, the phase diagram is found to be 
qualitatively different. For instance, some traffic states, 
which represent distinct phases in Ref. |^ff , make smooth 
crossovers to other traffic states without any sharp phase 
boundaries in between, implying that they are different 
limiting behaviors of a single dynamic phase. The in- 
vestigation is also performed for a large variety of initial 
conditions using two effective search methods. The con- 
ditions for the stable existence of the free flow, the RH 
state, and the traffic jam are examined. In some param- 
eter ranges, it is found that multiple dynamic phases can 
remain stable with respect to sufficiently small pertur- 
bations. In such parameter ranges, finite perturbations 
may induce transitions between those phases, resulting 
in metastability. Due to the presence of the on-ramp, 
the evolution process of the jam shows several different 
patterns and the phase boundaries for the formation of 
the jam are significantly modified. 

The paper is organized as follows. In the next section, 
we investigate the possible traffic phases for given val- 
ues of the upstream flux and the input flux through the 
on-ramp. Various new features are discovered, which are 
absent in homogeneous highways. We examine the con- 
ditions for the stabilities of the traffic jam and the RH 
state. The metastability among the free flow, the RH 
state, and the traffic jam is investigated, and the travel 
time distributions of the three states are compared. We 
also discuss the several different evolution processes of 
the jam due to the presence of the on-ramp. Based on 
the phase diagrams, we find that the on-ramp flux be- 
comes a more important factor for the formation of the 
traffic jam than the total flux, the sum of the on-ramp 
flux and the upstream flux. In Sec. Ill, we demonstrate 
analytically that our macroscopic model possesses non- 
trivial solutions which are indeed found in numerical sim- 
ulations. Finally, Sec. IV summarizes our results. 



II. PHASE DIAGRAMS OF TRAFFIC 
EQUATIONS WITH AN ON-RAMP 

In this work, we adopt the continuum model of the 
highway traffic flow proposed by Kerner and Konhauser 
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where p(x,t) is the local vehicle density and v(x,t) the 
local velocity. q- m (t)(p(x) is the source term representing 
the external flux through an on-ramp. The spatial dis- 
tribution of the external flux ip(x) is localized near x — 
(on-ramp position) and normalized so that qi n (t) denotes 
the total incoming flux. V(p) is the safe velocity that is 
achieved in the time-independent and homogeneous traf- 
fic flow. In Eq. (0), the second term on the right hand 
side represents an effective "pressure" gradient on vehi- 
cles due to the anticipation driving |f4| and the velocity 
fluctuations [p|Jf3[ , and the third term takes into account 
an intrinsic dampening effect that is required to fit the 
experimental data [p2| . Here t,cq, p are appropriate con- 
stants. The flux or flow, pv, is denoted below by either 
q or /. 

In order to investigate the effects of a single on- 
ramp, we use the open boundary condition. The up- 
stream boundary values of the density and velocity are 
fixed at p(x = —L/2,t) = p up and v(x = —L/2,t) — 
V(p U p), respectively. On the other hand, the values 
at the downstream boundary (x — L/2) are linearly 
extrapolated from their values at neighboring points, 
x = L/2 — Ax and L/2 — 2 Ax where Ax is spac- 
ing used in the discretization. The numerical simula- 
tions are performed using the two-step Lax-Wendroff 
scheme pj|. We choose the following parameters : r = 
0.5 min, p = 600 vehicles km/h, cq — 54 km/h, and 
V(p) = V Q (1 ~ p/p)/(l + E(p/p) 4 ) where the maxi- 
mum density p — 140 vehicles/km, Vb = 120 km/h, and 
E = 100 p3]. Concerning the discretization, spatial in- 
tervals of Ax = 37.8 m and time intervals of At = 10~ 4 
min are used. We choose the spatial distribution of 
the external flux as ip(x) = (2na 2 )~ 1 / 2 exp (— x 2 /2a 2 ) 
with a — 56.7 m. With this choice of parameters, crit- 
ical values arc pb = 21. 1 vehicles/km, p c i=25.3 vehi- 
cles/km, p C 2=62.3 vehicles/km, /b = phV{ph) = 2047 
vehicles/h, / c i = p c iV(p c i) = 2249 vehicles/h, and 
fc2 = Pc2V{Pc2) = 843 vehicles/h. The maximum flow 
that can be achieved in the time-independent homoge- 
neous flow is / max = max p {pV (p)} — 2336 vehicles/h. 
Below we are interested mainly in the low density regime 
(< pel), and thus we will use for brevity the subscript c 
in place of cl. 

In the real highway traffic, there are many kinds of 
noises which perturb the traffic out of its steady states. 
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The real traffic state is hence under an infinite sequence 
of perturbations and subsequent responses of dynamic 
states. Previous studies jl5|,[l8| on highway traffic with- 
out ramps however showed that many observed features 
of the traffic flow can be explained from the steady state 
properties of the continuum model without any noise. 
Motivated by previous successes, we will ignore noises in 
this paper. 

In the absence of noises, each dynamic phase of the 
traffic flow corresponds to a steady state, or equiva- 
lently an attractor of the nonlinear hydrodynamic model 
[Eqs. A steady state may exhibit complicated 

time dependences depending on the nature of the cor- 
responding attractor. We examine in this Section the 
linearly stable steady states (or phases of traffic flow) for 
a given upstream flux /(— L/2) = f up = p U pV(p up ) and 
the vehicle flux qi n (t) = /imp through the on-ramp at 
x = 0. Linearly stable states are, however, often unsta- 
ble to large perturbations and multistability can occur. 
Here it is worth emphasizing that the concept of the mul- 
tistability in dynamic systems is somewhat different from 
that in equilibrium systems. In equilibrium systems the 
free energy selects one particular state as a "true" stable 
state and other states become metastable. In dynamic 
systems, on the other hand, the free energy cannot be 
defined and the concept of the "true" stable state is not 
applicable. In this sense, all states are metastable and 
they all should be treated equally. 

Possible presence of multiple steady states makes it 
very difficult to search completely for all phases that are 
stable for given parameters, since it requires examina- 
tions of many different initial conditions. To search out 
all multiple steady states, we use two methods : One 
is to apply a triggering pulse to a steady state, for ex- 
ample, by changing the value of / rmp for a short time. 
For a sufficiently strong pulse, a transition to a different 
steady state can be induced, allowing the identification 
of a new steady state. The other is the adiabatic sweep- 
ing method. Starting with a given steady state for a 
particular set of the system parameters, one increases or 
decreases one parameter adiabatically. This way, one can 
find the range of the parameter values where a dynamic 
state remains stable. These two methods effectively sim- 
ulate a large variety of initial conditions. Using these, we 
investigate the steady states for given system parameters 
/ U p and / rmp . In particular, we concentrate on three rep- 
resentative values of / up , and for each of them construct 
a phase diagram for the entire range of / rm p- However, 
since too large a value of / rmp is unrealistic, we restrict 
our attention to the range / rmp < f££g = / max - / up . 

The values of / up studied in this work are chosen from 
the following considerations. In the previous studies of 
homogeneous highways without ramps, it was found that 
the flux /b provides an important boundary. Whereas 
the free flow is the only stable phase below /b, the traffic 
jam can be created above this value. In the presence of 
ramps, one can expect that the appearance of the traffic 
jam depends on whether the upstream flux / up is larger or 



smaller than /b. (Below we show that this expectation is 
not true, due to the nontrivial effect of an on-ramp.) This 
property motivated us to choose one representative value 
of / up in the range larger than /b and another smaller. 
We also choose a very small value for / up ( < C /&), which 
later reveals the importance of / rmp on the formation of 
the congested traffic. 

A. /up > f h 

The phase diagram of the traffic states for / up = 2119 
vehicles/h appears in Fig. 0(a). Here / r c mp is the criti- 
cal input flux through the on-ramp above which the free 
flow in the downstream of the on-ramp becomes linearly 
unstable. The critical on-ramp flux /j C mp is determined 
from 

/imp = fc — /up ? (3) 

where f c = f cl . For < / rmp < / c mp , the flux, both 
in the upstream and downstream, is lower than f c but 
higher than /b. Hence the traffic jam can be created 
from the free flow by triggering events but it does not 
appear spontaneously. The finite amplitude perturba- 
tion to / rmp generates a cluster, which grows to a traffic 
jam since the upstream flux is larger than /b. The traffic 
jam propagates to the upstream with its characteristic 
group velocity. 

When / rmp > / r c mp , the flux of the free flow in the 
downstream is larger than the critical flux f c and the 
free flow is linearly unstable with respect to long wave- 
length perturbations of infinitesimal amplitude. The 
growth of infinitesimal perturbations leads to sponta- 
neously formed clusters and as pointed out in Ref. Jl6| , 
complex sequences of traffic jams may appear in the 
downstream region. In a certain range of / rmp , we also 
observe that clusters form a periodic regular sequence. 
It turns out that this regular sequence is caused by the 
presence of the on-ramp, whose detailed discussion will 
be given in the next subsection. 

B- /up < /b 

In Fig. |l](b), we present the phase diagram for / up = 
1948 vehicles/h. The free flow can exist until / rmp reaches 
/rmp ; as in the previous subsection. When / rmp is smaller 
than 92 vehicles/h, the free flow (with a transition layer) 
is the only stable phase. 

For / rmp > 92 vehicles/h, we find another time- 
independent state beside the free flow, which is shown 
in Fig. ^(a). In our simulation, this new state can be 
generated from the free flow by applying the triggering 
pulse in / rmp for a short time. Far away from the on- 
ramp, the density and flow are homogeneous both in the 
upstream and downstream. Near the on-ramp, a local- 
ized cluster appears, which does not propagate in cither 



3 



direction but stays motionless. Due to this immobility, 
such state is named as the "standing localized cluster" 
(SLC) state in Ref. The immobility of the SLC state 
is in contrast to the situations without ramps where all 
inhomogeneities should propagate. Hence the property 
is due to a novel effect of the on-ramp. Another interest- 
ing property of the SLC state becomes manifest in the 
density-flow relations. Notice that the density-flow rela- 
tions [circles in Fig. ||(b)] measured at several locations 
near the on-ramp do not necessarily fall on the homoge- 
neous density-flow relation curve (p, pV(p)) [solid line in 
Fig. ||(b)] even though the relation at each measurement 
location remains stationary with time. More remarkably 
the circles lie in the linearly unstable density region. 

Incidentally, an experimental data which may be rele- 
vant to this has been reported ||. It was observed that 
when the traffic is in the stationary synchronized flow 
state, the density and flux can remain stationary during 
a relatively long time interval (2-5 min). Their station- 
ary values often lie in the linearly unstable density region 
and they form a two-dimensional area in the density-flow 
plane instead of falling on a single well-defined density- 
flow relation curve. In Ref. ||, the stationary values are 
interpreted as an indication of the spatially homogeneous 
traffic, and Helbing, Hennecke, and Treiber |H| proposed 
the HCT state as an origin of the stationary synchronized 
flow. The HCT state provides an explanation for the sta- 
bility of the traffic in the linearly unstable density region 
but it leads to the formation of the well-defined density- 
flow relation curve, failing to explain the absence of such 
a curve in the measurement. 

Present analysis of the SLC state raises an alterna- 
tive possibility. The SLC state shows that being station- 
ary does not necessarily imply the homogeneity, and it 
also explains the stability in the linearly unstable den- 
sity region. Furthermore it can explain the absence of 
the well-defined density-flow relation curve. We mention 
that upon the adiabatic variations of / up , / rmp and the 
external flux profile ip(x), the density- flow relation at a 
single measurement location can cover a two-dimensional 
area in the density-flow plane. These agreements raises 
an interesting possibility of an alternative explanation 
for the stationary synchronized traffic flow based on the 
SLC state. We judge however that it is yet premature to 
draw a definite conclusion from these agreements alone. 
Further experimental investigation of the stationary syn- 
chronized traffic flow is necessary. In the next Section, we 
demonstrate analytically that the traffic equations (|l]J|) 
do have the SLC state solution. 

As the on-ramp flux / rmp increases adiabatically, one 
finds the phase transition from the SLC state to the re- 
curring hump (RH) state [Fig. ^(a)] . In the RH state, a 
cluster, or a hump, does not remain stationary but moves 
back and forth in a localized region near the on-ramp. Its 
drift to far upstream is not allowed since the upstream 
vehicle density is lower than the boundary value pb- The 
RH state is investigated in detail in Ref. |l9| using the 
periodic boundary condition, and many interesting prop- 



erties are found such as the discontinuous transition from 
the free flow to the RH state induced by localized pertur- 
bations of finite amplitudes, hysteresis, gradual spatial 
transitions from the RH state to the free flow, and syn- 
chronized oscillations. These properties are identical to 
those of the synchronized flow (nonstationary type) [||,|| , 
and based on these common properties, the RH state is 
proposed as the origin of the synchronized flow. 

In addition to the properties of the RH state discussed 
in Ref. we investigate here the transition between 
the RH state and the SLC state. Our simulation shows 
that the transition from the SLC state to the RH state 
and the reverse transition occur at the same critical value 
°f /imp without hysteresis. We also examine the oscilla- 
tion amplitude of the RH state. The amplitude decreases 
to zero continuously as / rmp approaches the critical value 
[Fig. ||(b)]. Below the critical value, the hump does not 
oscillate and it becomes a standing cluster. These prop- 
erties suggest that these transitions are a result of the 
supercritical (or very weak subcritical) Hopf bifurcation 
p5[ of the SLC state to the RH state. 

These transitions between the SLC state and the RH 



state are not observed in the previous study |19|, where 
the adiabatic decrease of the ramp flux leads to the dis- 
continuous transition of the RH state to the free flow in- 
stead (Fig. 3 in Ref. (f{|). We attribute this difference to 
the different boundary condition adopted in this paper. 
Unlike the open boundary condition where / up and / rmp 
can be controlled independently, the periodic boundary 
condition used in (T^| is such that the increase (decrease) 
°f /imp is always accompanied by the decrease (increase) 
of / up since the average density of the total system is 
fixed. Therefore the "scanning" direction in Ref. |E| is 
different from that in this paper. 

We next discuss the traffic jam state. In homogeneous 
highways without ramps, the formation and propagation 
of the jam can not occur when the flux is smaller than 
/b. In the present case with an on-ramp, the flux / up 
in the upstream region is lower than /b while the flux in 
the downstream region can be controlled by / rm p- Thus a 
usual jam that consists of a single localized cluster should 
decay after they reach the upstream region. So in this 
sense, a usual traffic jam is not a steady state. Our inves- 
tigations show that a different type of traffic jams (Fig. ^) 
can occur even when / up < /b due to the nontrivial ef- 
fect of the on-ramp: Clusters are self-generated near the 
on-ramp repeatedly, forming a "train" of clusters mov- 
ing upstream. Although each constituting cluster decays 
during its upstream propagation, the train can still re- 
main stable provided that the decay rate is smaller than 
the self-generation rate, which is controlled by the extent 
of the inhomogeneity, / rmp , rather than by the upstream 
or downstream flux. Thus the stability limits of the new 
traffic jam do not coincide with those of the usual traffic 
jams [Fig. 0(b)]. Notice that this train structure is differ- 
ent from usual traffic jams in homogeneous highways. To 
indicate the structural difference, we will call this state 
the "oscillating congested traffic" (OCT) state. 
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We note that a traffic jam state very similar to the 
OCT state appears for / up > /b. Clusters are self- 
generated near the on-ramp repeatedly, forming a regular 
sequence of clusters. For / up > /b, however, each cluster 
does not decay in their upstream movement because of 
the high upstream density. 

The structure of the OCT state can be compared to 
the RH state. In both states, clusters appear recurrently 
near the on-ramp. In the OCT state, however, the area 
of the congested region expands with time, while in the 
RH state, clusters are localized. This difference is due to 
the larger size of clusters in the former. 

It is also worth mentioning that the structure of the 
OCT state shows an interesting crossover as / rmp varies. 
For small values of f lmp (close to the lower stability 
limit), the distance between the clusters is relatively large 
so that there exist homogeneous flow regions in between 
[Fig. ||(a)]. As /rmp increases, the distance between 
the clusters shrinks and for sufficiently large values of 
/rmp, the homogeneous regions between them disappear 
[Fig. ||(b)], and the clusters are "closely packed" inside 
the congested region. 

In Rcf. these structural differences are discov- 

ered using a different hydrodynamic model and the OCT 
states for small and large / rm p are identified as two dis- 
tinct phases. The former was called the "triggered stop 
and go" (TSG) flow and the latter OCT. In this paper, 
however, we find that these apparently different states 
transform smoothly to each other as / rm p is varied, with- 
out any signature of singularities. Thus we group these 
two states as a single dynamic phase in this paper. This 
difference between this paper and Rcf. J2l| may be due to 
the different models used, but presently we do not know 
the precise origin of the difference. 

We emphasize that in a certain range of / rm p, three 
phases, the free flow, the RH state, and the traffic jam 
(OCT) can coexist. In this metastable region of / rmp , 
small differences in the initial traffic condition may re- 
sult in quite different final states. We mention that in a 
recent measurement ||] , very similar initial states of the 
free flow are observed to undergo different phase tran- 
sitions either to the synchronized flow or to the traffic 
jam. Here we obtain the three phases from the traffic 
equations with a fixed parameter set. 

The difference between the three phases, the free flow, 
the RH state, and the traffic jam (OCT), is manifest in 
the travel time distributions which are shown in Fig. [|. 
In order to calculate the exact travel time distributions, 
we determine the trajectory of a vehicle, which is initially 
located at cc V eh(io), as follows: 

2Veh(i) = 2Veh(*o) + / dt' v(x vch (t'),t') . (4) 

From the trajectory of each vehicle, we obtain the vehicle 
travel time passing through the region from xq = —5 km 
to x% = 5 km. With the same / up = 1948 vehicles/h and 
/rmp = 222 vehicles/h, the travel time distributions of the 



three states show different behaviors. While it consists 
of a single peak for the free flow, those for the RH state 
and the traffic jam (OCT) show broad distributions due 
to the nonstationary nature of these phases. Also notice 
that in average, the travel time for the traffic jam (OCT) 
is greater than those for the free flow and the RH state. 



C /up « / b 

Fig. ^(c) shows the phase diagram for / up = 1497 ve- 
hicles/h (about 25% lower than /b = 2047 vehicles/h). 
The free flow remains linearly stable for / rm p < /rmp- 
In a narrow range of / rm p, 480 vehicles/h< / rmp < 492 
vehicles/h, the SLC state is found, and for / rm p > 492 ve- 
hicles/h, the OCT state is found. For this low upstream 
flux, however, the RH state does not appear. We find 
the critical value of / up below the RH state is absent is 
about 1872 vehicles/h. 

It is interesting to notice that the upper stability limit 
of the SLC state and the lower stability limit of the OCT 
state coincide within our numerical accuracy. We verified 
that upon the adiabatic increase of / rm p, the SLC state 
undergoes the phase transition to the OCT state, and 
upon the adiabatic decrease of / rm p, the reverse phase 
transition occurs, both at / rmp = 492 vehicles/h. This 
coincidence raises an interesting possibility of a close re- 
lation between the two phases. This possibility is also 
supported by the expansion rate of the congested region, 
which seems to approach zero smoothly as / rmp is re- 
duced to the lower stability limit of the OCT state. 

We next examine the evolution of the traffic jam state. 
Fig. H shows the evolution of the structure of the con- 
gested region as / rm p is increased. For relatively small 
/rmp, the structure is the same as in Fig. |]. As / rm p 
increases, however, a homogeneous flow region appears 
near the on-ramp, which expands with time [Fig. |6|(a)]. 
Hence the congested region is partitioned into an inho- 
mogeneous part and a homogeneous one. For an even 
higher value of / rmp (=730 vehicles/h), the inhomoge- 
neous part shrinks in length with time and after this 
transient process, the whole congested region consists of 
a homogeneous part [Fig. ^(b)]. In Ref. [^l), this state of 
traffic flow is named the "homogeneous congested traffic" 
(HCT) state and is identified as a distinct phase. 

Unlike Ref. plf , however, it is not so clear in our sim- 
ulations whether the OCT and HCT states are distinct 
phases. The distinction between the OCT and the HCT 
state is obscured further in our simulations by the pres- 
ence of an intermediate traffic state where both the OCT- 
likc inhomogeneous part and the HCT-like homogeneous 
part expand with time. We call this intermediate state 
the "mixed congested traffic" (MCT) state [Fig. |(a)]. 
As / rmp increases, the change from the OCT state to the 
MCT state and then to the HCT state seems to occur in 
a smooth way. We thus infer that the OCT, MCT, and 
HCT are different forms of a single jam phase. 
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We now focus on the HCT state. Notice that even 
though / U p < /b, the area of the congested traffic in- 
creases monotonically with the group velocity of the up- 
stream front about —7.68 km/h, which is considerably 
lower than the usual jam propagation velocity ~ —15 
km/h This monotonic widening of the congested re- 
gion is caused by the "blockage" effect of the ramp, which 
can be understood easily by recalling that in real high- 
ways, a large flux from the on-ramp can almost block the 
flow of vehicles on main highways. We also mention that 
the structure of the HCT state is identical to the traffic 
jam (shock) caused by a blockage [p6| , which implies that 
for large /rmp, the ramp works as a bottleneck. 

Interestingly, the density in the congested region lies 
in the linearly unstable region of the homogeneous flow, 
Pci(= Pc) < P < Pc2- Hence according to Refs. |l^ , |l5| , 
long wavelength fluctuations of even infinitesimal ampli- 
tude should grow in this region. In our simulations, we 
find that small inhomogeneities in the initial state indeed 
grow to form clusters. These clusters however disappear 
when they reach the upstream boundary of the congested 
region and the congested region becomes homogeneous 
afterwards. This result implies that inside the linearly 
unstable region, there exists a range of density where the 
homogeneous flow is convectively stable ^] (that is, the 
instability drifts away in one particular direction leaving 
the regions behind unaffected). The same observation is 
made in Ref. pi[ , where the HCT state is related to the 
stationary synchronized flow. As mentioned in the pre- 
ceding subsection, however, the density-flow relation in 
congested region of the HCT state lies on the single curve 
(p, pV(p)), which differs from the experimental observa- 
tions of scattered data [||. We mention that the HCT 
state may be very sensitive to the presence of noises since 
it is only convectively stable. 

It is instructive to compare the stability ranges in 
Fig. 0(c) with those in Fig. 0(b). The stability range of 
the SLC state lies entirely below f^ mp in Fig. 0(c) while 
it is mostly above f^ mp in Fig. 0(b). Also the stability 
range of the OCT reaches below f^ mp in Fig. 0(c) while it 
lies entirely above f^ mp in Fig. 0(b). This suggests that 
when the sum / up + / rmp is the same, the formation of 
the cluster is easier for larger / rmp and thus the actual to- 
tal flux level in the downstream is lower. Physically this 
tendency can be understood as resulting from the larger 
density gradient near the on-ramp when the relative por- 
tion of /rmp is larger. This trend is indeed observed in 
highway measurements |28| and is called the "capacity 
reduction" . 



III. ANALYTIC SOLUTIONS OF THE SLC AND 
THE HCT STATES 

In Sec. II, we showed that various forms of traffic flows 
occur near an on-ramp. In the case of the free flow and 
the usual traffic jam, they are affected by the on-ramp in 



minor ways and their properties are essentially the same 
as those without an on-ramp, which have already been 
investigated intensively [pl|. For other phases, however, 
the presence of the on-ramp is crucial and understanding 
of their properties are relatively poor. In this Section, we 
present analytic studies of two forms of traffic flow, the 
SLC state and the HCT state. 



A. Standing localized cluster (SLC) state 

The analytic examination of the SLC state is relatively 
simple since all time dependence disappears. This analy- 
sis also provides a good starting point for future analysis 
of the RH state and the OCT state since they are closely 
related to the SLC state as discussed in the preceding 
Section. Hence we present below the analysis of the SLC 
state in detail. 

In a homogeneous highway, inhomogeneities in the den- 
sity or the velocity always propagate. In the presence of 
an on-ramp, on the other hand, the numerical investiga- 
tion in the previous Section shows that inhomogeneities 
may form a standing cluster without propagation. Here 
we demonstrate analytically that the model [Eqs. (0J|)] 
indeed allows standing cluster solutions in the presence 
of an on-ramp. 

To obtain the SLC solution, one imposes 

m = m =0 - (5) 

By integrating Eq. (I) with respect to x, one obtains 



p(x)v{x) = /, 



rmp 



ip(x)dx + / U p = q(x) . (6) 



Since the function q{x) is completely determined for given 
/ rmp and /up, one can use this equation to express p(x) 
in terms of v(x). Using this, one can rewrite Eq. (0) as 
follows: 



dx 2 



Cg \ dv q(x) 



dx 



V 



q{x) 



v 

(7) 



where d/dx has been replaced by dj dx since all time de- 
pendence disappears. 

For further analysis, it is convenient to assume a partic- 
ular form of the influx profile ip{x). We take the localized 
influx limit and choose y>(x) — S(x). Then q(x) becomes 
/up for x < 0, /up + /rmp for x > 0, and q'(x) = f rmp 5{x). 
Then Eq. (0) can be decomposed into two separate prob- 
lems defined on two semi-infinite regions, x < and 
x > 0, with the matching conditions at x = , 



v(x)\ x=0+ 
dv 

x =o+ 



v(x)\ x=0 - 
dv 

dx *=o- 



pv(0) 



I rmp 



(8) 



G 



For each semi-infinite region, it is instructive to rewrite 
Eq. (0) as follows, 

p— = q s \l - — \ w V I — ) - v , 

ax \ w L J tv L V v / J 



where s=p(ositive) for x > and s=n(egative) for x < 0, 
and q p = f up + / rmp , q n = f up . Notice that after 
the variable transformations v — > y, x —> t, p — ► to, 
Eq. (|^) can be regarded as the equation of motion of 
a particle subject to a potential U s (y) where dU s /dy = 
(Qts/Ty)[V(q B /y) — y] and to the "strange" coordinate- 
dependent damping force. 

For the safe velocity V(p) adopted in this paper, 
Eq. (|9|) is highly nonlinear and it does not seem feasi- 
ble to write down solutions in a closed form. However, 
qualitative properties of the solutions can be still investi- 
gated by taking Eq. (^) as a set of flow equations defined 
on the phase space (v, w). 

Since the global structure of the flow is largely deter- 
mined from properties of fixed points, we first find fixed 
points of Eq. (||). Simple algebra shows that there are 
three fixed points, (v,w) — (0, 0), (v s i, 0), (v s2 , 0). The 
first one is unphysical since v = implies p — > oo. This 
unphysical fixed point appears since we set V(p) =0 for 
p > p, and we ignore this below. The other two come 
from the two solutions v s \, v S 2 (< v s i) of V(q s /v) = v. 
(It can be easily verified that for q s < / ma x, there are 
always two solutions, the larger one corresponding to the 
maximum point of the potential U s (y) and the smaller to 
the minimum point.) 

We are interested in the solutions where 

/ \ / «nl for x -> -oo , . 

V(X) — ► <^ r (10) 

I v p i tor x oo. v ' 

Thus (v n i,0) [point A in Fig. [j](a,b)] is the relevant fixed 
point for x < 0. By linearizing Eq. (S), it can be veri- 
fied that it is a saddle point. Then the flow [path 1 in 
Fig. 0(a,b)] that is associated with the unstable eigen- 
direction of the fixed point determines the entire flow in 
the semi-infinite region x < 0. Similarly (v pl ,0) [point 
C in Fig. |?](a,b)] is the relevant fixed point for x > 0, 
which is again a saddle point. The entire flow in the pos- 
itive semi-infinite region is then determined by the flow 
[path 2 in Fig. 0(a,b)] that is associated with the stable 
eigen-direction of the fixed point. 

In order to construct a legitimate solution from the 
path 1 and 2, one should join the two paths using the 
matching conditions [Eq. (g)]. It is convenient to regard 
the conditions [Eq. (g)] as a definition of a mapping de- 
fined in the phase space (v, w), from a point (v,w) to 
(v,w + Cp/rmp//- tu )- The effect of the mapping is shown 
in Fig. 0(a), where the path 1 is mapped to the path 3. 
The path 3 crosses the path 2 at a point denoted as F in 
the inset. Then one can construct a full solution by con- 
necting the curve AE with the curve FC. This solution 



exists for an arbitrary small value of / rmp and represents 
the free flow solution (with a transition layer [Q at the 
on-ramp) . 

Different solutions appear for sufficiently large / rmp . 
The mapping of the path 1 to the path 3 for a larger 
/ rmp is depicted in Fig. 0(b). Now the path 3 crosses the 
path 2 at three points F, H, I, which implies the pres- 
ence of three solutions. The solution associated with the 
crossing point F again corresponds to the transition layer 
solution. The two other solutions associated with the 
crossing points H and / correspond to the desired SLC 
solutions. Each solution provides the velocity profile for 
— oo < x < oo, from which the density profile can be 
obtained from p(x) = q s /v(x). Fig. |?](c) compares the 
density profile associated with the crossing point H with 
that from the direct numerical simulation of the traffic 
model [Eqs. with the Gaussian form of ip(x). No- 

tice that the two profiles are essentially identical except 
for a small difference in the ramp region, which arises 
due to the approximation of the input flux profile by a 
delta function. The density profile associated with the 
crossing point / can be obtained in a similar way. This 
solution, however, is not found in the direct numerical 
simulation, which suggests that this solution is linearly 
unstable. 

This analysis implies that the SLC solutions appear 
only for / rm p larger than a critical value. Precisely at 
the critical value, the two crossing points H and / coin- 
cide. For /rmp larger than the critical value, numerical 
simulations indicate that only one solution (one through 
the crossing point H) is stable and the other is not. Thus 
one finds that a turning point connecting the stable and 
unstable SLC solutions appear at the critical ramp flux. 

Above considerations use a specific form of an influx 
profile. We believe that the precise profile does not 
change the qualitative nature of above discussion. In 
fact, the value of the critical ramp flux for the SLC state 
obtained using the approximation (p(x) — S(x) is found 
identical within the numerical accuracy to that obtained 
from the direct numerical simulation of Eqs. (0,^|) (with 
tp(x) Gaussian). 



B. Homogeneous congested traffic (HCT) state 

When the on-ramp influx is added, we observe a new 
kind of the traffic jam: The downstream front is fixed at 
the on-ramp and the upstream front moves with a fixed 
group velocity. Between the downstream and the up- 
stream fronts, the congested region of the homogeneous 
flow is maintained [Fig. §|(b)]. We notice that the up- 
stream front is a steady structure in a proper reference 
frame. Below we show analytically that Eqs. |) pos- 
sess the HCT state solution. Since the congested region 
is homogeneous, one can split the discussion into two 
parts, one for the moving upstream front and the other 
for the fixed downstream front. For the upstream front, 
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the relevant equations of motion are 
d(pv) 



where 



dp 

at 

dv 
~di ' 



dx 
dv 
dx 







P 



V{p) 



dp 

dx 



P 



d 2 v 
dx 2 



(11) 
(12) 



Since this front is surrounded by wide regions of the ho- 
mogeneous flow both in the upstream and the down- 
stream, one can impose the following boundary condi- 
tions, 



p(x = 


— oo, i) 


= P- 


— Pup ) 


(13) 


v(x = 


— oo, i) 


= v_ 


= V{p-) , 


(14) 


p(x = 


+00, t) 


= P+ 




(15) 


v[x = 


+oo, t) 


= v + 


= V(p+) . 


(16) 



Here p- is the density of the far upstream region and 
p + of the congested region Also the spatial coordinate 
is chosen in such a way that x = +oo corresponds to a 
location deep in the congested region instead of the far 
downstream in the original equations p|,2[. 

Let us assume that Eqs. ( pd[ ) and (|12) allow a steady 
state solution that satisfies the boundary conditions 
Eqs. ([|), (|l|), (|l|), and (|l|). Here the steady 
state means that in a proper reference frame, all time- 
dependence disappears. We perform the change of the 
reference frame: 



x 
t' 



x 
t , 



Vgt 



(17) 
(18) 



and neglect all time-dependence in this new reference 
frame to find 



q 



d4_ 

dx 1 
, dv 
dx 1 



= 



P- 



V(p) 



• dp 
] dx' 



ft 



d 2 v 
dx 12 



(19) 
(20) 



where q' = pv — pv g . We can determine the constants 
q' and v g from the boundary conditions Eqs. jl^), (14), 
(|l5|), and (filf). From the condition that the in- flux to 
and the out-flux from the front should be the same in 
the primed reference frame, we obtain 



q' = 



P+v+ ~ P-V- 
P+ ~ P- 

P+P-( v - - " 



P- 



Using Eq. (pi), one has 



P = 



<1 



(21) 
(22) 

(23) 



Plugging in this expression into Eq. (EOh, one finds 



d 2 v 
'dx' 2 



_(V - Vg)'' 



- 1 



dv 



— +F(v } q\v g ) = 0, (24) 



F(v,q',v g ) 



T(V - Vg) 



V 



(25) 



Eq. (^4j) is again an equation of motion of a particle 
subject to a conservative force —F and the coordinate- 
dependent damping force. The boundary conditions 
Eqs. (|l3|), (|l4| ) ensure that F(v±; q', v g )~0 automatically. 
Since u_ > v+, the root v~(v+) corresponds to the po- 
tential maximum (minimum). Thus the nature of the 
stationary state is clear in the particle motion analogy. 
At x = — oo, the particle is at the unstable maximum 
point v — V—. As "time" x increases, the particle slides 
down the hill and after some time, it settles down at 
v = V+ due to the friction, provided the damping coeffi- 
cient in Eq. (Q) remains positive. This particle motion 
describes the upstream front of the HCT state [Fig. ||(b)]. 

Before we begin the next analysis of the downstream 
front of the HCT state [Fig. |(b)], one remark is in or- 
der. In the direct numerical simulation of Eqs. [0, 0, the 
vehicle velocity in the congested region of the HCT state 
is fixed for given / up and / rm p- However, in the above 
analysis of the upstream front, where / up is used to fix 
V- = v up , the value of v + is still a free parameter. We 
show below that this degree of freedom should be used 
to allow the downstream front solution. 

Analysis of the stationary downstream front is very 
similar to the analysis of the SLC state in the preceding 
subsection. Using the condition of the stationarity and 
integrating Eq. (]]]), one obtains Eqs. (|5|) and (^), respec- 
tively. Adopting the approximation ip(x) = S(x), one 
recovers the matching conditions at the ramp [Eq. (^)] 
and the set of the flow equations [Eq. ([)[)] . The value of 
q n should be to allow a continuous connection of 

v(x) to the upstream front solution and q p = q n + / rm p- 
The asymptotic behavior of v(x) far away from the on- 
ramp should be chosen differently as follows: 



v(x) 



Vn2 
Vpl 



for x 
for x 



■ — oo 
oo , 



(26) 



where v p i and v n 2 are defined in the same way as in the 
preceding subsection, and x — > — oo corresponds to the 
region deep in the congested region. 

Notice that for x ~ * — oo, v(x) approaches v n 2 instead 
of v n i since v+ corresponds to the smaller of the two so- 
lutions of V(q n /v) — v [or since the corresponding den- 
sity p+ lies on the descending slope of the homogeneous 
density- flow relation]. This difference in the asymptotic 
behavior results in a qualitative change. It can be ver- 
ified through the linearization of Eq. (|9|) that (w n 2,0) is 
a stable fixed point. Then we find v(x) — v n 2 for the 
entire semi-infinite region x < 0, which should be con- 
trasted to the preceding subsection. Near the fixed point 
(upl, 0), on the other hand, the situation is similar to the 
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preceding subsection and the flow in the positive semi- 
infinite region becomes a continuous path, like the path 
2 in Fig. |f(a,b). 

To construct a full solution of the HCT downstream 
front, the separate solutions for x < and x > should 
be joined using the matching condition [Eq. (Q)]. Since 
the velocity for x < is constant, the matching condi- 
tions reduce to 



v(x) 


\x=0+ 


= v n2 , 


dv 




r 2 


dx 


x=0+ 





This matching conditions can be satisfied only when 
the phase space trajectory for x > passes the point 
{v n 2,<$/(jJ>v D 2)f Tmp ). For given v n2 and / rmp , the path 
in general does not pass the point except for a particular 
value of v a 2 — v+- This tuning thus fixes the free param- 
eter as mentioned before. In general, u+ is a function 
of the influx profile tp(x) since the matching conditions 
depend on ip(x). 

Fig. H shows the density profile (solid line) of the HCT 
state obtained from this matching method for <p(x) — 
S(x). It is essentially identical to the profile (dotted line) 
obtained from the numerical simulation of Eqs. (0,^|), ex- 
cept for the larger density peak at the on-ramp caused 
by the approximation tp(x) = S(x). We confirm from the 
simulation result that the group velocity of the moving 
front is consistent with —7.68 km/h which is given by 
Eq. ( pl| ) and the damping coefficient q'[c%/(v — v g ) 2 — 1] 
is positive in the congested region. 

IV. CONCLUSION 

The traffic equation with a source term representing 
the on-ramp influx of a highway displays a variety of 
novel traffic flow states not present in the homogeneous 
equations. To understand the role of the source term, 
we map out the phase diagram using the continuum traf- 
fic model [Eqs. proposed by Kerner and Konhauser 
. In our numerical simulation, we use the open bound- 
ary condition, which allows one to handle the single on- 
ramp without using very large system sizes. Due to pos- 
sible presence of multiple metastable states, detailed sim- 
ulation is carried out for a limited number of represen- 
tative values of the upstream flux / up and for the whole 
range of the on-ramp flux / rmp . Various traffic states are 
identified and characterized. The phase diagrams thus 
obtained are summarized in Fig. |l|. It is found that an 
inhomogeneous but stationary traffic state (SLC) can ap- 
pear near the on-ramp due to / rm p- This state is related 
to the recent measurement of the homogeneous synchro- 
nized flow. The on-ramp also generates a new kind of 
traffic jam, which can appear even below the stability 
limit of the usual traffic jam in homogeneous highways. 
The structure of the new traffic jam varies qualitatively 
with /rmp- The capacity reduction due to / rmp is also 



observed. In a certain range of / mp , the free flow, the 
RH state, and the new traffic jam can all coexist so that 
the free flow can undergo phase transitions either to the 
RH state or to the traffic jam state. Analytic investiga- 
tions are also performed and two nontrivial solutions of 
Eqs. (0J|) are found. These solutions describe the SLC 
state and the HCT state. 
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FIG. 1. The phase diagrams for / up > / b (a), / up < / b 
(b), and / up < f h (c). Here f h = 2047 vehicles/h, f c = 2249 
vehicles/h, and / max =2336 vehicles/h while f^ mp = fb — /up, 
/rmp = fc - /up, and = /max - /up. (a) The phase dia- 

gram for /up = 2119 vehicles/h where /r mp = 130 vehicles/h, 
and /^J^=217 vehicles/h. For / rmp > /r mp , the traffic jam 
state is generated spontaneously, (b) The phase diagram for 
/up = 1948 vehicles/h where / r b mp = 99 vehicles/h, / r c mp = 301 
vehicles/h, and f™p = 389 vehicles/h. The metastable re- 
gion among the free flow, the RH state, and the OCT state 
extends from / rmp = 206 vehicles/h to 238 vehicles/h. (c) The 
phase diagram for / up = 1497 vehicles/h where /!m P =550 ve- 
hicles/h, / r c mp = 752 vehicles/h, and / r ™ a p x = 839 vehicles/h. 
As the on-ramp flux increases, the traffic jams with different 
structures appear. The lower stability limit of the HCT state 
is denned as the value of / rmp above which the inhomoge- 
neous part disappears. In a similar way, one may define the 
upper stability limit of the OCT as the value of / rmp below 
which the homogeneous part does not expand with time. No- 
tice that between these two stability limits, there exists an 
intermediate range of / rm p for which both the OCT-like part 
(inhomogeneous part) and the HCT-like part (homogeneous 
part) expand with time. This intermediate state is called the 
"mixed congested traffic" (MCT) state. The stability limit 
between the OCT and MCT state is / rmp = 603 vehicles/h 
and that between the MCT and HCT state is / rmp = 730 
vehicles/h. 
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FIG. 3. (a) The evolution of the RH state for / up = 1948 
vehicles/h and / rm p = 222 vehicles/h. The hump moves back 
and forth in a localized region near the on-ramp. (b) The 
evolution of the RH state for / up = 1948 vehicles/h and 
/rm P =130 vehicles/h. Near the lower stability limit of the 
RH state, which coincides with the upper stability limit of 
the SLC state, the oscillation amplitude of the RH state is 
very small. 
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FIG. 2. (a) The spatial density profile of the SLC state 
for /up = 1948 vehicles/h and / lmp = 121 vehicles/h. The 
on-ramp is at km. The profile does not change with time, 
(b) Circles: The density-flow relation for several positions 
near the on-ramp. That denoted as U(D) is the data in the up- 
stream (downstream) homogeneous region. Each data point 
is stationary with time, and upon adiabatic variations of / up 
and the external flux profile tp(x), it covers a two-dimensional 
area. Solid line: The curve q = pV(p). 
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FIG. 4. (a) The spatiotemporal evolution of the density of 



the OCT state for / up = 1948 vehicles/h and /„ 



222 vehi- 



cles/h. The values of / up and /rmp are the same with those in 
Fig. |(a). The RH state in Fig. |(a) and the OCT state in this 
figure are independent metastable states, (b) The spatiotem- 
poral evolution of the OCT state for / up = 1948 vehicles/h 
and /rmp = 381 vehicles/h. The increase of / rm p generates 
the "closely packed" clusters. Similarly to (a), each cluster 
decays at far upstream from the on-ramp. 
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FIG. 6. The evolution of the MCT (a) and HCT state 
(b) for / up =1497 vehicles/h. (a) / rmp =635 vehicles/h. The 
congested region consists of the homogeneous part and the 
inhomogeneous part, (b) / rmp =794 vehicles/h. The inhomo- 
geneous part is not present. The upstream front moves with 
a fixed group velocity. 



FIG. 5. The travel time distributions of the free flow (solid 
line), the RH state (dashed line), and the traffic jam (OCT) 
(dotted line) for / up = 1948 vehicles/h and / rmp = 222 vehi- 
cles/h. The travel time distributions are obtained by following 
10 5 trajectories of vehicles through the region from -5 km to 
5 km. 




12 





sz 

0) 

> 



110 



90 



70 



50 



30 



10 



I 


I I 

\ 


A J 




i i i 



-5 km 



km 
x (km) 



5 km 



(c) 

FIG. 7. The flow in the phase space (v, w) (a) with / rmp 
slightly below the lower stability limit of the SLC state and 
(b) with /rmp slightly above the lower stability limit. The 
path 1 represents the flow that is associated with the un- 
stable eigen-direction of the fixed point (« n i,0) (point A), 
and the path 2 the flow that is associated with the stable 
eigen-direction of the fixed point (v rmp i , 0) (point C) . The 
path 3 represents the image of the path 1 through the map- 
ping ([s]) . The insets are the enlargements of the dotted regions 
at right, (c) The density profile (solid line) of the SLC state 
obtained analytically for / up =1948 vehicles/h and / rmp =121 
vehicles/h. The marks A, G, H, C represent corresponding 
points in Fig. |^(b) . The density profile obtained from the nu- 
merical simulation of Eqs. g|) is also given for comparison 
(dotted line, vertically shifted 40 vehicles/km). Notice that 
the density jump at the on-ramp (x = 0) in the solid line 
(due to the approximation tp(x) = S(x)) is smoothed out in 
the dotted line. 
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FIG. 8. The density profile of the HCT state obtained an- 
alytically (solid line) under the approximation tp(x) = 5(x) for 
/up = 1497 vehicles/h and /rm P =762 vehicles/h. The density 
profile obtained from the numerical simulation of Eqs. (pl,0) 
is also given for comparison (dotted line, vertically shifted 50 
vehicles/km). Small differences between the two are due to 
the different choices of tp{x). 
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